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Abstract 

We show that near a focus-focus point in a Liouville integrable 
Hamiltonian system with two degrees of freedom lines of locally con- 
stant rotation number in the image of the energy- momentum map are 
spirals determined by the eigenvalue of the equilibrium. From this rep- 
resentation of the rotation number we derive that the twist condition 
for the isoenergetic KAM condition vanishes on a curve in the image 
of the energy-momentum map that is transversal to the line of con- 
stant energy. In contrast to this we also show that the frequency map 
is non-degenerate for every point in a neighborhood of a focus-focus 
point. 
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1 Introduction 



Consider a neighborhood of an equilibrium point of a Hamiltonian system 
with two degrees of freedom. In local symplectic coordinates (x,y,p x ,p y ) £ 
R with standard symplectic structure £1 = dx A dp x + dy A dp y the Hamilto- 
nian is H(x,y,p x ,p y ), where the equilibrium is taken to be the origin. The 
Hamiltonian vector field in these coordinates is 

& H = JdH, J=(\ J) (1) 

or in components 

._8H._dH.dH.dH 

X dp x ' ^ dp y x dx' ^ y dy 

Since an equilibrium of a dynamical system is a zero of the vector field one 
has dH(0) = 0. The critical value H (0) of the critical point is assumed to 
be 0. The equilibrium is characterized by the eigenvalues A of its lineariza- 
tion Jd 2 -fT(0). If A = a + iuj is an eigenvalue, so is A, —A, and —A. Real pairs 
of non-zero eigenvalues ±a are called hyperbolic, non-zero pure imaginary 
pairs ztio; are called elliptic, and quadruples of eigenvalues with non-zero 
real and non-zero imaginary part ±a ± iu; are called loxodromic. Near a 
loxodromic equilibrium there exist symplectic coordinates (x,y,p x ,p y ) such 
that the quadratic part H2 of H is (231 EE| 

H 2 = aJi + ujJ 2 , Ji=xp x + yp y , J 2 = xp y - yp x . (3) 

Introducing z = x + iy and p z = p x — \p y gives Q = $t(dz A dp z ), and 
Ji — 1J2 = zp z . Outside z = this extends to a multivalued symplec- 
tic coordinate system (In z,zp z ), and the flow of H 2 is easily found to be 
(z,p z ) = (exp(At)zO) exp(— Xt)p z o). The quadratic Hamiltonian has a two- 
dimensional stable eigenspace spanned by the eigenvectors of eigenvalues 
with negative real part, and a two-dimensional unstable eigenspace spanned 
by the eigenvectors of eigenvalues with positive real part. In each eigenspace 
the dynamics is that of a focus point in the plane. The invariant manifolds of 
the equilibrium point of the full Hamiltonian are tangent to these eigenspaces 
at the equilibrium. In general the stable and unstable invariant manifolds 
will have transverse intersections, and the system is non-integrable. 

Here we are concerned with the dynamics near an equilibrium with lox- 
odromic eigenvalues of a Liouville integrable Hamiltonian system with two 
degrees of freedom. Under some additional hypothesis (see below) we call 
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such an equilibrium focus-focus point. Let L be an independent second inte- 
gral that commutes with H, i.e. the Poisson bracket {H, L} = 0,(JdH, JdL) 
vanishes. The system is called Liouville integrable if the energy-momentum 
map F = (H, L) is regular almost everywhere, and {H, L} = 0. The values 
in the image of F are denoted by c = (h,l). The Liouville- Arnold theo- 
rem states that near any compact connected component of regular points 
T c of F _1 (c) there are so called action-angle coordinates (9%, 62, Ii, I2) with 
Q = d6i A dl\ + d#2 A d/2 and I\ and I2 are commuting constants of motion 
with periodic flows on (and near) the two-dimensional torus T c C F _1 (c). 
In these coordinates the Hamiltonian is a function of 1% and I2 alone, and 
the quasiperiodic flow of H is the solution of 9\ = dH/dli = cui(I), i = 1, 2. 
In other words, the Hamiltonian vectorfield 3Cn is a linear combination of 
the periodic flows of the actions 

^H = ^ h +U 2 ^I 2 . (4) 

The coefficients of the linear combination are the frequencies oJi{I) of the 
angles 9i, i = 1,2. The frequencies uj = (cji,^) in general depend on the 
actions I = (Ii,^)- The map from actions I to frequencies co(I) is called 
the frequency map. 

The KAM theorem (see e.g. p]E3) asserts that under small perturba- 
tions the invariant torus T c with frequencies uj(I) persists if 1) the frequency 
map is non-degenerate and 2) the frequencies u>(I) are Diophantine. In the 
isoenergetic KAM theorem persistence of invariant tori is considered for 
the same energy h. Then a torus is characterized by its rotation number 
W(I) which is the frequency ratio \u\ : 102] € MP 1 . The iso-energetic non- 
degeneracy condition dW/dl 7^ 0, where the derivative is taken at constant 
energy, replaces the non-degeneracy condition on the frequency map, see 
e.g. pp. In a semi-global Poincare section transversal to 3>h on T c the 
iso-energetic non-degeneracy condition ensures that the Poincare map is a 
twist map near the invariant curve corresponding to T c . The twist condition 
dW/dl 7^ implies that the rotation number changes between neighboring 
tori. An invariant torus T c for which the twist vanishes is called a twistless 
torus, for short. To check the non-degeneracy of the frequency map in an 
example requires tools from complex analysis, see e.g. Horozov's work on 
the spherical pendulum JOJ- General results in the neighborhood of critical 
points of the energy momentum map are due to Knorrer jl.Sj and Nguyen 
Tien Zung |24| . We are not aware of general results about the iso-energetic 
non-degeneracy condition; the spherical pendulum has again been treated 
by Horozov |llj . 
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The twist condition is only a sufficient condition for persistence, and 
KAM theorems with weaker conditions exist, see ^7j. However, the per- 
turbed dynamics is quite unusual when the unperturbed twistless torus is 
resonant, i.e. when it has a rational rotation number. A resonant torus with 
twist breaks into a Poincare-Birkhoff island chain, see e.g. |151 1181 12]. A 
twistless resonant torus instead breaks into two island chains, and near the 
collision of these chains interesting dynamics occurs on so called meander- 
ing invariant curves, see [I'll |31 1191 IS] and the references therein. Our goal 
is to show that such dynamics occurs near a loxodromic equilibrium of a 
generically perturbed integrable system. More precisely we will prove the 
following 

Theorem 1. In every integrable Hamiltonian system with two degrees of 
freedom and a focus-focus singularity with loxodromic eigenvalues there ex- 
ists a regular torus with vanishing twist for each value of the energy close to 
the critical one. For this energy all other tori close enough to the singular 
fibre have non-vanishing twist. 

This theorem might look surprising when compared to a result by Nguyen 
Tien Zung |24j : He showed that near a focus- focus point the Kolmogorov 
non-degeneracy condition is always satisfied, i.e. on almost all tori the 
frequency map is non-degenerate. Zung did not provide a way to find the tori 
for which the non-degeneracy condition fails; but, as we show in section 0J 
our techniques improve his theorem by showing that the frequency map is 
non-degenerate for all regular tori close enough to the singular fiber. 

It is well known that the Kolmogorov condition and the Twist condi- 
tion are independent, and our results show that whenever there is a focus- 
focus point the twist-condition is violated while the frequency map is non- 
degenerate. 

2 Rotation Number 

Let the critical value of the equilibrium be (h, I) = (0, 0). For simplicity, we 
assume that the corresponding singular fibre -F -1 (0) contains only a single 
critical point m. We then need to assume that the component of -F _1 (0) 
containing m is compact. Such a singularity is called of focus-focus type. 
Since we are interested in tori close to the singular component of i ? " 1 (0) 
we may disregard other connected components. Then (0, 0) is an isolated 
critical value in the image of the energy-momentum map F restricted to 
those tori. Additional structure appears near the loxodromic equilibrium by 
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the assumption that the system be integrable. Eliasson showed that the 
momenta J\ and J2 of the quadratic normal form are the components of a 
momentum map J = (Ji, J 2) of the full system near the focus-focus point in 
some symplectic coordinates. This means that we may assume that both H 
and L (and not only their quadratic parts) are functions of J\ and J2 near 
m: 

H = $ o J and L = * o J, . (5) 

where (7 = ($, \P) is a local diffeomorphism of M? near the origin. Recall that 
(J2J) is valid only near the focus-focus point m. However this gives a way of 
extending the momentum map J to a whole neighborhood of the singular 
fiber simply by enforcing J = L). 

Since H is loxodromic the coefficient a in Q is non-zero hence <9i$(0) 7^ 
0. Then by the implicit function theorem the additional integral L can be 
chosen equal to J2. 

Observe that J2 has a 2ir periodic flow near m, and this still holds in a 
neighborhood of the singular fiber, as was shown in Therefore J2 is a 
generator of an S 1 action on the neighboring tori and hence is equal to one 
action, say I2. 

Instead of using the energy-momentum map F it is advantageous to use 
the momentum map J instead, and whenever the Hamiltonian enters the 
calculation to use 

In |22| it is shown using only the momentum map J without the Hamil- 
tonian (or, equivalently, with Hamiltonian equal to J\) that the 2n periodic 
flow SCi x of the first action is given by 

27^=^^+72^ (6) 

where the periods 73, i = 1, 2 satisfy 

n(j) =a x (j)- M(ln(), 
r 2 (i)=a 2 (i) + 3f(lnC), 

and o"j(j), i = 1,2, are smooth and single- valued functions near the origin. 
Here the point j = (ji, j'2) in the image of the momentum map J is identified 
with the complex number C = Ji + ij2- The smooth contribution comes 
from the dynamics far away from the focus-focus point, while the singular 
contribution is obtained from the flows of J\ and J2 near the focus-focus 
point, see for the details. The points (71,72) and (0,2rr) form a basis of 
the period lattice of the foliation. Using this result we obtain 
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Lemma 2. The rotation number W(j) near a focus-focus point is 

2nW(j) = -A(j)to(\n () - 3(ln C) + a{j) (8) 
where A(j) is smooth and determined from the Hamiltonian H = o J by 

^) = ff(A m = <») 

and a (J) is a smooth (and single-valued) function near the origin. 
Proof. The Hamiltonian vector field 3Cn is obtained from JHJ) as 

x H = d^&j, + a 2 $jr j2 . (io) 

Now eliminating Xj x with the aid of © and using J 2 = I 2 gives 

^ = — JT/, + (d 2 $ - -di$) 3£ l2 . (11) 

Comparing with (JIJ the frequencies tUj can be read off, and the rotation 
number as a function of j is 

w(j) = = - tsCj)) • (12) 

Finally using (JJJ) gives the result with a = Ao\ — <7 2 ■ From (j3J) we see that 
A(p)=u/a. □ 

The fact that 3:(ln£)/27r is multivalued and increases by 1 upon com- 
pleting a cycle £ = exp(i0) is a manifestation of the fact that a simple 
focus-focus point has Monodromy with index 1, see [5| I14[ l25 | 12 lj . 

The functions t% that determines the period lattice depend on the folia- 
tion alone, while A{j) is determined by the Hamiltonian alone. In order to 
understand the behaviour of W(j) near the origin we now prove the following 

Theorem 3. There is a local diffeomorphism of M 2 near the origin which 
is C 1 at the origin and C°° elsewhere which maps the level sets of W to 
the integral curves of the equation £ = —AC- In particular these curves are 
spirals determined by the eigenvalue A of the focus-focus point. 

Proof. We use the notation of lemma 2. The strategy is to change the 
coordinates such that the level sets of W in the new coordinates satisfy a 
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linear differential equation. Let ip be the change of variables defined in polar 
coordinates ( = ji + U2 = pe ld by 

/ p = /O)M(0)exp(-a(j)M(0)) \ 

\e=e )■ 

Writing A(j) = A(0) + (j,A(j)), where A is smooth, we have ip(j) = 
f(j)j, where / : M 2 — > R locally near the origin is given by 

f(j)= p ^Mj))/A(0) eM _ a(j)/m y 

We see that / is continuous at the origin, with /(0) = exp(— a(0)/A(0)), 
and that the partial derivatives of / are of order In p at the origin. Hence (p 
is C 1 at the origin, with Jacobian 

dp(0) = /(0)ld = exp(-a(0)/A(0))Id. 

Under this diffeomorphism, the rotation number is simplified and we 
only need to study the level sets of 

2nW o p-\j) = (-A(O)K(InC) - S(lnC)). 

Forget now the tildes. Viewing If as a local Hamiltonian for the standard 
canonical structure dj'2 A dji , we compute the levels sets of W as the integral 
curves of the associated Hamiltonian vector field: 4-(ji,j2) = {(hW, —d\W). 
Rescaling the time s by 27r|C| 2 a and rewriting the dynamical system in 
complex form gives the equation £ = —XQ. 

From Lemma 2 we know that W(j) diverges when j approaches the 
origin. These statements are compatible because different branches of the 
complex logarithm are involved. After each complete turn of the spiral the 
rotation number jumps by one. Alternatively one could view the rotation 
number as globally defined on the Riemann surface of the complex logarithm. 

□ 

Recall that as coordinates in E? we may use j = (ji, j'2) or c = (h,l), 
and that they are related by I = ji and h = Therefore the spirals in 

the image of the momentum map J will be mapped to spirals in the image 
of the momentum map F. Near the origin the transformation between the 
two is the linear map (h,l) = {j\a + 

3 Vanishing Twist 

In the isoenergetic problem fixing the energy H = h gives a smooth curve 
^/i in the image of the momentum map. The twist condition is obtained 
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from W(j) by differentiating along this curve. Since 9i$(0) = q / we 
may apply the implicit function theorem and the curve ^ is the graph of a 
function ji(j 2 )- The slope of this curve near the origin is given by implicit 
differentiation of $(ji,j 2 ) = h which yields dj\/dj 2 = — ^2$/9i$ = —A(j). 
Therefore the twist is found to be 



= diW$± + d 2 TU = -A(j)diW + d 2 W . (13) 



9h 

Expanding the remaining derivatives gives 

2irT(j) = -A 2 d 1 r 1 + 2Ad 2 Ti - d 2 r 2 - r x Ad x A + nd 2 A (14) 

where d\T 2 = d 2 T\ has been used, which follows from the fact that both peri- 
ods can be obtained by differentiating a single action function, see |22| . The 
main contribution near the origin comes from differentiating the complex 
logarithm in (|7j). which produces terms proportional to l/|j[ 2 . Therefore we 
now define T(j) = 27r|j| 2 T(j) and have the following 

Lemma 4. T is a C 1 map at the origin that satisfies 

t(0)=0, dit(O) = A 2 (0)-l, d 2 f(0) = -2A(0) (15) 

Proof. Using where o~i is smooth we find that the last two terms in (|14|). 
TiAd\A and T±d 2 A, are both of the form f(j) +g(j) In \ j\ with smooth / and 
g, while the remaining first three terms are of the form /(j) + g(j)/\j\ 2 for 
some smooth / and g. Since |j| 2 ln|jj is of class C 1 at the origin the first 
statement follows. A simple computation gives 

t(j) = A 2 h - 2Aj 2 - h + 0(|i| 2 In |j|) , (16) 

and the result follows. □ 
Since the derivatives of T cannot both vanish at the origin, T~ 1 (0) is a 
C 1 curve through the origin. By Lemma 21 the equation for the tangent of 
this curve at the origin is 

(A(0) 2 - l)h = 2A(0)j 2 (17) 

Using h = j\ct + j 2 oj at the origin to express this tangent in (h,j 2 ) gives 

J 1 + a 2 , 

h = uj— ^j 2 - (18) 

When u> ^ the curve T(0) therefore transversally intersects the curves 
{h = const} near the origin, and we have proven Theorem 1. Note that the 
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condition oj ^ is only needed in the last step. Therefore also for cu = 
twistless tori exist near the origin, but transversality is lost. Therefore 
twistless tori might exist only for positive or only for negative values of h, 
but possibly not for both. 

Different choices of actions are possible. The semi-global S 1 action L = 
J2 = h is unique up to sign e = ±1, but the other action I\ can be changed, 
I\ = eli + kl2, k € Z, so that new and old actions are related by a uni- 
modular transformation. The frequencies uo change into uj\ = eu>\ and u)2 = 
eu>2 — kuji, so that the rotation number transforms as W = W—ke. Therefore 
a different choice of actions amounts to a different choice of sheet of the 
Riemann surface of the complex logarithm. This argument also shows that 
the vanishing twist is independent of the choice of actions. 

When there are more than one critical point in the singular leave similar 
results can be derived, based on modified formulas (JJJ) whose derivation is 
sketched in [22]. 

As an application of the above results one can consider the integrable 
normal form of the Hamiltonian Hopf bifurcation, see |2L)l |5| and the refer- 
ences therein. In [Jj the rotation number in the compact case was computed 
explicitly in terms of elliptic integrals. Expansion of these integrals at the 
focus-focus point gave the above results (JHJ) and (|18j) for the first time in this 
special case. Pictures of the spirals of constant W can also be found in [7]. 
But the method employed could not deal with the case in which the normal 
form has a non-compact singular leave. And even with a compact singular 
fibre in a more complicated system this approach might lead to hyperelliptic 
integrals, and their expansion at the singular point would be quite difficult. 
A prominent example of an integrable system with a focus-focus point is the 
spherical pendulum, see e.g. [3]. In this case the eigenvalues are degenerate 
because oj = 0, and the spiral of Theorem |21 degenerates into a star. 

4 Kolmogorov Condition 

Near a focus-focus point we can improve Zung's theorem |24j . 

Theorem 5. On regular tori close to the focus-focus singular fiber, the 
Jacobian determinant of the frequency map admits the following asymptotic 
expansion (recall that t\ is of order In 




(19) 
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In particular this shows that the Kolmogorov condition is uniformly satisfied 
for every torus near the singular fibre of a simple focus-focus point. 



Proof. First notice that the Jacobian of the frequency map does not depend 
on the choice of action variables and hence does not care about monodromy; 
in our case, this means that it does not depend on the determination of the 
complex logarithm in (JJJ) • Thus the leading order of this Jacobian can easily 
be calculated. We need to compute 

doj dui dJ 
~dl = ~dJ~dI' 

From @ we know that by definition 
oJ\ 0J2 

Since J2 = I2 we easily find that det dl/dJ = t\/2it. The frequencies can 
be read off from (|11|) : they are 

UJl = $1, 0J 2 = $2 $1- 

Tl Tl 

The subscripts of * denote partial derivatives. Using that t\ is of order 
In I j I and that the first derivatives of n have all leading order l/\j\, one 
easily computes: 

d x u} 1 = -2^*^ + 0(1^1) =°(rrr) 

Tl Wl 

3 2 uj x = -2vr*i^ + 0(l/ri) = 0(r7T72) 



' J l-?1 T l (Oil) 

a (T1T21 - T 2 Tii) 1 V > 

dito 2 = 2 -$i + 0(l =0-— 

Ti \j\n 

r, (TlT22 - T2T12) , . . 1 . 

5 2 cj 2 = 2 $i + 0(l) = 0(— — ) 

Ti bin 

Here the second index of r denotes a partial derivative. At leading order t\ 
can be replaced by In \j \ . The leading order of the determinant det duj / dJ 
is obtained by taking the determinant of the matrix in which only the terms 
that are given explicitly in (j20j) are kept, provided one shows a posteriori 



that the result has the expected order of jj^pr = rjp ln a j^i ■ The result is 
9a; < I >2 1 

det in = MniT22 - ti 2 2 )^- + 0( 



aj --V.xx.- '12, r 3 ' -V |i|r 2. 
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Using the form of T\ and T2 we find 



det 7— = Tiir 2 2 - t\ 2 = ---2 + 0(l/\j\) 



so that 





2 



2 



1 



-2vr 



+ 0{ 



Ijpln 4 



Since by hypothesis the leading term has indeed the required order. 

Returning to the true actions / the final result is 



Remark. In the process of finishing this paper we became aware of a preprint 
by Rink |16| in which the determinant of the Jacobian of the frequency map 
was calculated in a very similar way. 
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